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Abstract 

In this article, we derive source integrals for multipole moments in 
axially symmetric and static spacetimes. The multipole moments can 
be read off the asymptotics of the metric close to spatial infinity in a 
hypersurface, which is orthogonal to the timelike Killing vector. Whereas 
for the evaluation of the source integrals the geometry needs to be known 
in a compact region of this hypersurface, which encloses all source, i.e. 
matter as well as singularities. The source integrals can be written either 
as volume integrals over such a region or in quasi-local form as integrals 
over the surface of that region. 

1 Introduction 

In general relativity, there were several definitions of multipoles proposed. Since 
this theory is non-linear, it is, however, by no means obvious that this is at 
all possible. Thus, it is not surprising that in early works multipoles were 
only defined in approximations to general relativity that lead to linear field 
equations and allow a classical treatment. The most definitions in this direction 
and beyond were covered in Thome's review [Tj. 

From the 1960s on, new definitions in the full theory of isolated bodies^ 
started emerging. These definitions of multipole moments can roughly be di- 
vided into two classes. In the first, the metric (or quantities derived from it) 
are expanded at spacelike or null-like infinity. We will call these asymptotic 
definitions or asymptotic multipole moments. Amongst these are the definitions 
of Bondi, Metzner, Sachs and van der Burgh (BMSB) [2][3], Geroch and Hansen 
(GH) 4,5 , Simon and Beh£] (SB) Janis, Newman and Unti (JNU) [7]|8], 
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1 This means that all sources (matter and black holes) are located in a sphere of finite 
radius and the spacetime is assumed to be asymptotically flat. A precise meaning is given in 
Se.pOl 

This approach reproduces the GH multipole moments. 



Thor ne [T] , the ADM approach [9] and the Komar integrals 10 , for reviews 
sec []][llf There scope of applicability varies greatly. Whereas the GH multi- 
pole moments arc defined only in stationary spacetimes the BMSB, JNU, Thorne 
and ADM definitions hold in a more general setting. The Komar expressions 
for the mass and the angular momentum on the other hand require stationarity 
and stationarity and axially symmetry, respectively. Higher order multipoles are 
not defined in the Komar approach. Despite their conceptual differences, Giirsel 



showed in 12 the equivalence of the GH and Thome's multipole moments in 
case the requirements of both definitions are met. Additionally, the mass and 
the angular momentum in the GH, Thorne, ADM and Komar approach can be 
shown to agree. 

In the second class fall multipoles that are determined by the metric in a 
compact region. Dixon's definition in [13] falls in this class. These multipoles 
are given in the form of source integrals. However, it is not yet known how 
they are related with the asymptotic multipole moments. A main application 
of these multipoles is in the theory of the motion of test bodies with inter- 
nal structure. But for test bodies it is obvious that there cannot be any such 
relation between Dixon's multipole moments and the asymptotic multipole mo- 
ments. Furthermore, Dixon's definition is in general not applicable if caustics 
of geodesies appear inside the source, i.e., if the gravitational field is too strong 



compared to a characteristic radius of the source. Ashtekar et al. defined in 14 



multipole moments of isolated horizons. These are also source integrals and 
require only the knowledge of the interior geometry of the horizon. In 14 , it 
was also shown that the so defined multipoles of the Kerr black hole deviate 
from the GH multipoles. This effect becomes more pronounced the greater the 



rotation parameter. In both cases 13 14 , it is interesting to find the relation of 



these multipole moments to asymptotically defined GH multipole moments and 
their physical interpretation. Source integrals might proof useful in this respect. 
Other definitions, in particular for the quasi-local mass and the quasi-local 



angular momentum, can be found in 15 . Here we aim at source integrals for 



all multipole moments in static spacetimes for arbitrary sources, i.e., we want 
to express asymptotic multipole moments by surface or volume integrals, where 
the surface envelopes or the volume covers all the sources of the gravitational 
field. That such source integrals can be found, is not at all obvious. This is 
because of the non-linear nature of the Einstein equations yielding a gravita- 
tional field, which acts again as a source. To overcome the principle difficulties, 
we will focus here on static and axially symmetric spacetimes. In this case, the 
vacuum Einstein equations can be cast in an essentially linear form. Addition- 



ally, they allow the introduction of a linear system 16-18 . The latter might 
seem superfluous, if the the former holds. However, we will derive our source 
integrals relying solely on the existence of the linear system and the applica- 
bility of the inverse scattering technique. This will allow us in future work to 
apply the same formalism to stationary and axially symmetric isolated systems. 
This is especially relevant for the description of relativistic stars, where many 
applications are. A generalization to spacetimes with an Einstein-Maxwell field 
exterior to the sources seems also feasible. In all these cases, multipole moments 
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of the GH type exist. 

Source integrals will prove useful in many respects, in particular in the search 
of global solutions of Einstein equations describing figures of equilibrium or 
relativistic stars, for recent efforts see, e.g., 19-22 and references therein. For 
example, an exterior solution can be constructed for a known interior solution by 
employing the source integrals to calculate its multipole moments. From these 
the exterior solution can be completely determined. This yields an exterior 
solution, which does not necessarily match to the interior solution. However, if 
it does not, this construction shows that there is no asymptotically flat solution, 
which can be matched to the given interior. Conversely, possible matter sources 
can be analyzed for given exterior solutions and their multipole moments. The 
latter approach is also of astrophysical interest, since often only the asymptotics 
of the gravitational field and the asymptotic multipole moments are accessible 
to experiments. Source integrals can be employed to, e.g., restrict the equation 
of state of a rotating perfect fluid from observed multipole moments. 

Furthermore, source integrals can be used to compare numerical solutions, 
analytical solutions and analytical approximations by calculating their multi- 
pole moment^J see for example 24 . This can also be used to approximate 



the vacuum exterior of a given numerical solution by an analytical one, which 
exhibits the correct multipole moments up to a prescribed order, see 25 . The 
merit of source integrals lies in the fact that they determine the multipole mo- 
ments using the matter region only, which captures the internal structure of the 
relativistic object and is usually determined with high accuracy. Additionally, 
source integrals provide the means to test the accuracy of numerical methods, 



which are used to determine relativistic stars, cf. 26 , by calculating the mul- 
tipole moments in two independent ways: Firstly, using the asymptotics and, 
secondly, using the source integrals. This will give also a physical interpretation 
to possible deviations. 

The paper is organized as follows: In Sec. [5J we introduce the different 
concepts used later, i.e., the GH and the Weyl multipole moments as well as the 
inverse scattering technique. Sec. [3] is devoted to the derivation of the source 
integrals and includes the main results, Eq. ((23} -([24}. In Sec. [3) we will discuss 
some properties of the obtained source integrals. 



2 Preliminaries 

In this section, we will repeat the notions that are needed in the present paper. 
Note, that we use geometric units, in which G = c = 1, where c is the velocity 
of light and G Newton's gravitational constant. The metric has the signature 
(— 1, 1, 1, 1). Greek indices run from to 3, lower-case Latin indices run from 1 
to 3 and upper-case Latin indices from 1 to 2. 

3 For difficulties in extracting the multipole moments from a given numerical metric, see, 
e.g, [23]. 
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2.1 The line element and the field equation 

We consider static and axially symmetric spacetimes admitting a timelike Killing 
vector £ a and a spacelike Killing vector rj a , which commutes with £ a , has closed 
timelike curves and vanishes at the axis of rotation. If the orbits of the so 
dchncd isometry group admit orthogonal 2-surfaces, which is the case in vacuum, 



in static perfect fluids or static electromagnetic fields, see, e.g., 27 , then the 
metric can be written in the Weyl form: 

ds 2 = e 2k - 2U (dp 2 + d( 2 ) + W 2 e~ 2U d^ 2 - e 2U dt 2 , (1) 

where the functions U, k and W depend on p and £0 Note that the metric 
functions U and W can be expressed by the Killing vectors: 

e 2C/ = -£ Q r, W 2 = -W&C". (2) 

The Einstein equations can be inferred from the ones given in 129 . Since we 
will not specify the matter, we give only a complete set of combinations of the 
non-vanishing components of the Ricci tensor: 



A(2)c/ + w {u > pW ' p + UxWx) = e2fc ~ 4t/i? «< 

W, pp - W xc + 2 (k x W x - k p W, p + W (U 2 p - U 2 C )) = W (i? cc - R pp ) , 

W x k, p + W, p k, ( ~ 2WU X U P - W :PC = WR pC , (3a) 

WA^W = e 2k - iU W 2 R u - e 2k R vv , 

-2A^k = (R pp + R CC ) - e 2k - iU R tt - ^R, 



where A^ n ^ = {j^ji + + ~§c^j ■ The fourth equation implies that we can 

introduce canonical Weyl coordinates (p, £) with W = p via a conformal trans- 
formation in vacuum or in all domains, where A' 2 ^M^ = holds, including, e.g., 
dust. After this coordinate system is chosen, we drop the tilde again. The 
remaining coordinate freedom is a shift of the origin along the symmetry axis, 



which is characterized by p = 0. Eqs. (3a I simplify in vacuum to the well-known 
equations 

A (3) t/ = 0, 

k <c = 2pU, p U k p = p((U iP ) 2 ~(U, c ) 2 ), 



where A^™) is defined analogously to but with canonical Weyl coordinates. 
The last two equations determine k via a line integration once U is known. 



This k automatically satisfies the last equation in (3a I. Hence, only a Laplace 



4 That this is not the general case, can be seen from the conformally flat, static and axially 
symmetric solutions in [28] . The field equations are not imposed there, but they can be used 
to define a (rather unphysical) stress energy tensor. 
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equation for U remains to be solved. Therefore, the Newtonian theory and 
general relativity can be treated on the same formal footing. We will do so here 
and highlight the difference in Sec. |4j The disadvantage of using the canonical 
Weyl coordinates is that they cannot necessarily be introduced in the interior 
of the matter, where we have to use other (non-canonical) Weyl coordinates. 




Figure 1: The different surfaces S and Si, volumes V and Vi are depicted for 
a certain matter distributions (ellipsoid and torus). The curves B and A are 
relevant in Sec. |2.3| Not that this serves only as illustration and that under 
a certain energy condition and separation condition static n-body solutions do 
not exist, for details see 30 



The physical setting, which we want to investigate, is as follows: We have 
in a compact region V of a hypersurface given by t = const, several sources, 
cf. Fig. [T] The surface of V will be denoted by S = dV and is assumed to be 
a topological 2-sphere. Note that inside V it will in general be not possible to 
introduce canonical Weyl coordinates. The sources can be matter distributions 
supported in Vi with Si — dVi- If there are black hole sources with horizons 
%i, we assume that there is vacuum at least in a neighborhood of Hi, which 



necessary for static black holes, see, e.g., 31 . Hence, we can define a closed 
surfaces Sf* in the vacuum neighborhood, which encloses only the black hole 
with T-Li^\ For relativistic stars, we can choose V = Vi, in which case we assume 
that Israel's junction conditions are satisfied across S — Si, i.e., that there are 
no surface distributions. In fact we will assume this for all surfaces Si and S^ . 
The only other restriction on the matter is that we want to be able to introduce 
Weyl's line element, see ([I]). 



5 Similarly, we can introduce Si around other singularities, if they have a vacuum neigh- 
borhood. 
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2.2 Geroch's multipole moments 

Isolated gravitating systems are described by an asymptotically flat spacetime. 
The precise meaning of this is defined below and used to define Geroch's mul- 
tipole moments. Let us denote our static spacetime by (M, g) with the metric 
g a p and by V a hypersurface orthogonal to £ Q endowed with the induced met- 
ric h a p = —tpi-ygaf) + £a£,3, for which we will use subsequently Latin indices. 
For the definition of Geroch's multipole moments, asymptotic flatness of V is 
sufficient, see [I]: 

Definition 1 V is asymptotically flat iff there exists a point A, a manifold V 
and a conformal factor fl 6 C 2 {V), such that 

1. V = V\JK 

2- h a f, = fl 2 h a i, is a smooth metric ofV 

3. £1 = D a Q — and D a Dt,Q — 2h a (, at A, where D a is the covariant deriva- 
tive in V associated with the metric h a b- 

Let us define the potential t/j 

S 2 2 

— ~ (3) 

which is also a scalar in V. If we introduce the Ricci tensor W ab ' built from the 
metric h a b then the tensors P ai ...a n can be defined recursively: 

P = $ 



Pa 1 ...a n — C 



DP (n-l)(2n-3) jya 



(5) 



where C[A ai ...a n ] denotes the symmetric and trace-free part of A ai ... 0n . The 
tensors P ai ...a n evaluated at A define Geroch's multipole moments: 

M 01 ... 0n =P 01 ... 0B | A . (6) 

The degree of freedom in the choice of the conformal factor reflects the choice 
of an origin, with respect to which the multipole moments are taken, see |2J. 

Up to now, only stationarity was required. If we choose axially symmetry 
the multipole structure simplifies to 



m r = —P a 
r! 



(7) 



where the z a is the unit vector pointing in the direction of the symmetry axis 
and the scalars m n define the multipole moments completely. Hence, we will 
refer to them as multipole moments, as well. 
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Fodor et al. demonstrated in 32 that the multipole moments can be 



obtained directly from the Ernst potential on the axis or in the here con- 
sidered static case from U at the axis. If we expand U along the axis, i.e., 

oo 

U(p = 0,C) = X) ^ |CI~ r ~ 1 ! we characterise the solution by the set of con- 

r=l 

stants U {r \ The result in [32] relates now the m r with the U^ r \ i.e., if the 
are known, we can in principle obtain the m r . Thus, we will limit ourselves to 
discussing mainly the lf( r \ Although the relation m r (U^') can be obtained in 
principal to any ord er, up to now only the mo, . . . , rnio were explicitly expressed 



using 



m 



the U (r \ 



SCO 



32 



We will give here only the first four for illustration: 



mi 



m 2 



3 



;y(3) 
(8) 



Eq. ^ shows that the mass dipole moment, E/W, can be transformed away, 
if the mass, U^°\ is not vanishing. For a general discussion, further references 
and expressions of the center of mass, see [33] . 

In 34 a method to obtain the m r was proposed, which could help to over- 
come the non explicit structure of m r (t/w). Also the pure 2 r pole solutions 



in 35 could proof useful in this respect. 



2.3 The linear problem of the Laplace equation 

Lastly, we shortly review the linear problem associated with the Laplace equa- 
tion. Although the equations involved are fairly simple, we decided to use this 
technique, because it is readily generalizable to the stationary case. 

In the more general case of stationarity and axially symmetry, the Einstein 



equation, i.e., the Ernst equation, admits a linear problem, see, e.g., |16}|18| and 
for a recent account 36 . In the static case[^] the linear problem reads 



il + X)(:,o-. a,: .!•')•' (!)) 



where z = p + iC,, the spectral parameter A = J x + \ z z , K £ C and a bar denotes 
complex conjugation. The function a depends on z, z and A. The integrability 
condition of Eqs. ([9| is the first equation in Eq. (3b). 



Next we will repeat some known properties of without proof. For details 



we refer the reader to 36 . There are four curves of particular interest A ± , B 
and C, cf. Fig. [T] The axis of symmetry is divided by V in an upper and lower 
part A + and A~ , respectively. The curve B generates S by an rotation around 
the axis and is given by a restriction of S to <p = 0. Thus, we will subsequently 
refer to A and B as curves in a p, C-plane. Lastly, C describes a half circle with 
sufficiently large radius connecting A + with A~ . 



°The formulas can easily inferred from 36] by setting gtip = 0. 
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Along A ± and C Eq. @ can be integrated. This yields for a suitable choice 
of the constant of integration 



(10) 



a(X = 


+1,P = 


o,C) 


= F(K)e 2U{p 


<r(A = 


-l.P = 


o,C) 


= 1, 


a(X = 


+1,P = 


o,0 


_ e 2C/(p=0,C) ; 


a(X = 


-hp = 


o,0 


= F(if). 



The function F : C -> C is given for K € M with (p = 0, Q = K) e A ± by 



F(K) = r 2U{P= ° X = K) We** 
1 ' i^(P=o.C=ic) ( ,F)e^-- 



The integration along 23 is the crucial part for our considerations in the next 
section. 



3 Source integrals Geroch's multipole moments 



Let us assume that the line element is w ritten in canonical Weyl coordinates 

The scalars U and W, cf. 



2.1 



are 



in the exterior region V c , cf. Sec 

also scalars in the projection V. Furthermore, U and W are supposed to be 
continuously differentiable in the vacuum region including S and, thus, B. Then 
we can consider the linear problem ([9| also along B (after a projection to the 
ip = plane): 



U A s A + 



i + AW + i (i 



A ) U A n A 



(12) 



where (s A ) = { s p,s^) = (§f , f s ) and (n A ) = (n",n c ) = (-^,gf) denote the 
tangential vector and the outward pointing normal vector to the curve B : 
s e [sjv,ss] — > (p( s ),C( s )), respectively. The parameter values spj/s give the 
"north/south" pole Pn/Si i- e -, [p = 0, C = (n/s)> ci - Fig- PI Note that the 

ahz 



tangential and the normal vectors are not necessarily normalized allowing an 
arbitrary parametrization of B. 

Eq. ( 12 1 constitutes an ordinary differential equation of first order with the 
boundary conditions as given in ( 10 ) assuming (0, K) £ A + U A~ . As such it is 
an over determined system, which corresponds to the integrability of Eqs. ([9]). 



However, Eq. ( 12 ) is readily integrated and the compatibility condition of the 



boundary conditions reads 



1 



U(0,K) =-(U(0,( N ) -U(0,(s))+ 

i 7 (13) 
- / ((A" 1 + A) U, A s A + i (A" 1 - A) U, A n A ) ds. 
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Eq. ( 13 1 determines the axis values of U from the Dirichlet data and the Neu- 
mann data along an arbitrary curve B, which is sufficient to obtain the entire 
solution U in V c . 

The multipole moments follow from an expansion of Eq. ( 13 ) with respect 
to K~ l . Let us denote by the expansion coefficient to order |X| _r_1 of a 

oo 

function f(K), which is constant at infinity i.e., f(K) = f^ r '\K |~' 1 . The 

r=-l 

coefficients 7V| r) = (A -1 + A)M and N { I ] = i(A _1 - A)M depend still on (p,() 



and satisfy the equations 



(r) 



N 



(r) 



N 

r(r) 



N^l = 0, 

r(r) 



P 



(14) 



Eqs. ( 14 1 follow directly from the form of the spectral parameter A, cf. after 
Eq. ([9]). This expansion is only valid for p 2 + £ 2 < 00. Furthermore, N± and 
their radial derivatives evaluate at the axis to 

r(r) 



N y _>(p = 0,() =0 Vr>-1, 
7V|- 1) (P = 0,C) = 2, 
ivj r) (p = 0,C) =0 Vr>0, 

^(P = o,C)=o 

iV^(p = 0,C) = -2C r Vr>0, 



(15) 



The zeroth order of Eq. (13) implies together with Eq. (15 1 that 1 ' = 0. 
This is also required by asymptotic flatness, which is assumed in the derivation 
of (10). Solving the Eqs. (14) yield after a lengthy calculation N± for r > 
everywhere: 



L§J 



*< r) = £ 

fc=0 



2(-l) fe+1 r!p 2fe+1 C 



r-2k 



A k {k\) 2 {r -2k)\ 



(16) 



k=0 



-2k-l 



2(-l) k+1 r\p 2k+2 C 
4 fc (fc!) 2 (r-2/c-l)!(2/s + 2) 



IV ' • 

Note that 

p 



is well-behaved also for p — > for r > 0. 



For orders r > Eq. (|13|) yields together with Eq. ( 16 ) the following line 

pole moments 

J r ' I (n^U a s a + N { _ r) U A n A ) d 7 , (17) 



integrals defining Weyl's multipole moments 

1 



713 
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where s A and h A are the normalized vectors s A and n A , respectively, and d7 
denotes the proper distance along the path 70, which runs along B from the 

(r) 

north to the south pole. The functions N± and U are to be read as functions 
of ( P (s),C(s)). 



Eqs. ( 17 ) are already expressions of the kind we are searching for, since they 
determine the multipole moments from the metric given in a compact region, 
i.e., they are quasi- local. But we also can rewrite these multipole moments as 
volume integrals justifying the term source integrals even better. The main 



obstacle for doing so is that Weyl's multipole moments are given in Eq. (17 1 
using canonical Weyl coordinates. Hence, neither the coordinate invariance of 
these expressions is transparent nor is obvious how to continue p and C to V. 
However, U can be expressed by the norm of the timelike Killing vector, cf. Eq. 
([2]), which can easily be continued to the interior. 
Let us introduce the 1-form 

Z a = e^sW'PW-^t 6 , (18) 

where e Q/ 3 7 5 denotes the volume form of the static spacetime. Z a is closed in V c 
and hypersurface orthogonal in the entire spacetimes. Thus, we can introduce a 
scalar potential Z with Z^ a = XZ a , where the scalar X equals 1 in the vacuum 
region. In canonical Weyl coordinates in V , the potential Z has the trivial 
form Z = C + Co- Because we did not fix the origin of our Weyl coordinates, 
e.g., the value of Cat, we can set the constant Co = without loss of generality. 
This integration constant is exactly the freedom we need to change the origin 
with respect to which the multipole moments are defined allowing us to change 
to the center of mass frame. W = p and Z = C in V c and they are defined 
everywhere. Hence, we can use these two scalars to continue p and C hito V. 
Note that W iQ , and Z a are orthogonal everywhere and have the same norm in 



the V c . The line integral (17) in the covariant form reads now 



Uir) = \f K r \w( S ),Z( S ))U a s a + N { _ r) (W(s),Z( S ))U, a n a d 7 . (19) 



The dependence on W(s) and Z(s) will be suppressed in the following expres- 
sions. 

Along B the functions W and Z satisfy 

W, s = Z <n , W, n = -Z,„, (20) 

which is a consequence of the field equation and the choice of canonical Weyl 
coordinates. After an partial integration, we can rewrite the line integral as 



surface integral using the axial symmetry and Eq. ( 20 1 , which yields 

Uir) = Ljw ( N - )U ' A - N +W Z * U + N +, Z W t nU) AS. (21) 
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We denote by d<S the proper surface element of S, Si or Sj 1 , respectively. Since 
we ruled out surface distributions, Israel's junction conditions imply that Eqs. 
( 21 ) can be understood as integrals over the 2-surface S as seen from the exterior 



or the interior, see 37 . 

Using Stoke's theorem and Eqs. ^ we obtain (see Fig. [T] and the end of 



Sec. 



2.1 



Jj{r) = 



for a description of the Sf 1 and Vj) 



1 

8^ 



N^iW Z) f a f^ 



W 



N[%(W,Z)U 



z>° 
W 



N 



(r) 



u 



U r Z (w,z) w (w> a w, a - z< a z <a ) 



dV- 



-Y 



II" 



(22) 

Note that the normal vector points outward at Sf 1 and dV is the proper 
volume clement of V or V;, respectively. The covariant derivative with respect 
to the metric g a p is denoted by a semicolon. The field equations in the vacuum 
region (3b I imply that the integrand vanishes there. Hence, we can write Weyl's 
multipole moments as the contributions of the individual sources to the total 
Weyl moment 



N ( ^(W,Z) 



W 



N { + ^ W (W,Z)U 



N 



(r) 



u 



+ ,wz(W,Z) w (W a W, a 



Z' a Z„) 



-E 

8tt ^ 



W 



u. 



N { 1] w ZmU- 



dV- 



Z^ a 
W 



N^ Z W, A U)<\S 



(23) 



(r)H 



The integrals in (23) are the source integrals or quasi-local expressions for the 
asymptotically defined Weyl moments and, thus, for the asymptotically defined 
Geroch-Hansen multipole moments. Note that the second derivatives of W in 



Eq. (24) can be expressed by the energy momentum tensor using Eq. (3a 



How the contributions of the black holes are related to the definitions of the 
multipole moments of isolated horizons in 14 will be investigated in a future 
work as well as the relation of the source integrals of Vj to those of Dixon given 
in, e.g., pi. 
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Since the transformation from Weyl's multipole moments to Geroch's is non- 
linear except for the mass and the mass dipole, cf. Q, there is no linear 
superposition of the multipole contributions of the individual sources to the 
total Geroch multipole moments. Hence a mixing of the contributions of the 
individual sources takes place. 

Of course, Eq. (23) can again be rewritten as a sum of surface integrals: 



8tt ^ 



w 



N^yWnU ) dS+ 



8tt 



s n W \ 



N%Z A U + N^ z W A U)dS. 



(24) 



It is obvious that our choice of continuation of p, £ into V is not unique and 



affects greatly the form of Eqs. (21 )- ( 24 ) , though not the value. In fact, any 



C 1 extension of the scalars W, Z from V c could be chosen. Thus, depending 
on the applications, other choices might be more appropriate. 



4 Properties of the source integrals 



In this section, we will discuss the consequences of Eq. (17) in more detail in 
the Newtonian and the general relativistic case. In the Newtonian case, Eq. 
(17) comprises the well known multipole definitions as we will show in the next 



section. 



4.1 The Newtonian case 

Suppose U is a solution to the Laplace equation in V c (cf. Fig 1), then it follows 
by virtue of Green's theorem from the Dirichlet and Neumann data: 



U{x) = - 



1 

4tt 



G(x,y) 



dU(y) 
dy a 



-U(y) 



dG(x,y) 
dy 



v 



(25) 



where G(x,y) denotes an arbitrary Green's function, dS y a surface element of 
the boundary S and x G V c . The integration is over y £ S and n a is the inward 
pointing with respect to V c unit normal to S at y. 
Eq. 

This 



(25) is equivalent to Eq. (21) in case one restricts x to the axis and 
makes an expansion in for the special choice G(x,y) = — ■ 

yields surface integrals in the form of (21) with the factor e u set to 1, a flat 
space surface element and with W = p and Z — Q. Thus, the multipole moments 



still contain all the information as Eq. (25) with the difference that the latter is 



not accessible for stationary, axially symmetric and isolated sources in general 
relativity. This is the reason, why we chose the approach using the linear system. 
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If Stoke's theorem is applied for these surface integrals, we arrive at 

[/M = i- / —A^UdV = J / — ^dV, (26) 
St: J p 2 J p 

v v 

where we made use of the Poisson equation for the Newtonian gravitational 
potential, A^U = 47r/z, with a mass density p. These are, of course, the usual 
multipole moments in source integral form up to a sign. This justifies the term 



'source integrals' also for the equivalently obtained expressions (23) in curved 
spacetime. A comparison with the well-known formulas of Newtonian theory 
shows that 

N ( _ k) = -2pr k P k {cos9), Vr > (27) 

with polar coordinates (r, 9) defined as usual: p — r cos 9, Q = r sin 9. P k 
denote the Legendre polynomials of the first kind. In the general relativistic 
case, these N± ' depend on the extensions of p, ( into the interior and thus 
become polynomials in the scalars W and Z. 

4.2 The general relativistic case 

In this section, we will discuss some of the properties of the quasi-local volume 



integral given in Eq. (22). The Einstein equations are non-linear and contain 



already the equations of motion (Bianchi identity). Hence, we could not expect 



a result like ( 26 ) , which depends only on the mass density. We rather find 



source integrals (23) containing terms that are not expressed explicitly by the 
matter distribution (all but the first term) . However, all of these terms vanish 
in vacuum. They also vanish in matter distributions, for which we can choose 
W = p in V. In those case the source integrals have the same form as in the 
Newtonian case. 

Of course, the first multipole moment coincides with the Geroch mass and, 
hence, must coincide with the well-known Komar mass. With Eq. (l8|) and Eq. 



we have N { " ] = o and ivl 0) = -2W such that only the first term of the 



integrand in Eq. (|22j) remains: 

M =hVl ^Mr™ + i^L eUuMS 

This is, of course, exactly Komar's integral of the mass in static spacetimes. 
The black hole contributions can also be cast in the standard form: 

M s n = ~ f eVU^dS = ^ / e aPl sC' P - (28) 

Although the main goal of this paper is to present the derivation and defini- 
tion of the source integrals, we will give here a short application. We show that 
static, axially symmetric and isolated dust configurations do not exist. This is 
an old resullQ but can easily be recovered using source integrals. This demon- 

7 For more general non-existence results for dust, see also 381440 . 
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strates also how these quasi-local expressions can be employed. Static and axi- 
ally symmetric dust configurations are characterized by the energy-momentum 
tensor in Weyl coordinates 

T ab = ^S'Jl (29) 
The Bianchi identity implies U a — in V and, thus, at S. This yields together 



with the quasi-local surface integrals for the Weyl moments Eq. ( 19 1 



U {r) = 0. (30) 

Thus, the system has no mass or any other multipole moment, which implies 
flat space in the vacuum region. This is clearly a contradiction to a dust source 
with positive mass density. 



5 Conclusions 

We have derived in this article source integrals or quasi-local expressions for 
Weyl's multipole moments and, thus, for Geroch's multipole moments for axially 
symmetric and stationary sources. These source integrals can either be written 
as surface integrals or volume integrals. A priori, one could not expect to 
find any kind of source integrals at all, because of the non-linear nature of 
the Einstein equations. That this is possible in the here considered setting, 
seems not to be due to the staticity and axially symmetry and the peculiarly 
simple form of the field equations. But rather a linear system must be available 
offering a notion of integrability of the Einstein equation. Thus, it appears 
feasible to find source integrals not only for stationary and axially symmetric 
isolated systems, which describe vacuum, but also electrovacuum close to spatial 
infinity. These generalizations will be investigated in future work. 

It should also be clarified, how the source integrals are connected to the 



already known source integrals for isolated horizons 14 . In 14 it was shown 
that the source integrals characterize the horizon uniquely. However, they do 
not reproduce the GH multipole moments of a Kerr black hole. In our approach, 
the agreement of the source integrals and the asymptotically defined Weyl or 
Geroch multipole moments is given by construction. Therefore, these source 
integrals might prove useful for identifying the contributions to the multipole 
moments, which yield the discrepancies between the isolated horizon multipole 
moments and those of Geroch and Hansen. 
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